Saddles, Twists, and Curls: Shape Transitions in Freestanding Nanoribbons 
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Efforts to modulate the electronic properties of atomically thin crystalline nanoribbons requires 
precise control over their morphology. Here, we perform atomistic simulations on freestanding 
graphene nanoribbons (GNRs) to first identify the minimal shapes, and then employ a core-edge 
framework based on classical plate theory to quantify the width dependence in more general systems. 
The elastic edge-edge interactions force ultra- narrow ribbons to be flat, which then bifurcate to 
twisted and bent shapes at critical widths that vary inversely with edge stress. Compressive edge 
stresses results in twisted and saddle shapes that are energetically indistinguishable in the vicinity 
of the bifurcation. Increasing widths favor the saddle shapes with (longitudinal) ribbon curvatures 
that vary non-linearly with width and edge stress. Positive edge stresses result in a flat-to- curled 
transition with similar scalings. At large widths with negligible edge-edge interactions, rippling 
instabilities set in, i.e. edge ripples and midline dimples for compressive and tensile edge stresses. 
Our results highlight the utility of the core-edge framework in developing a unified understanding of 
the interplay between geometry and mechanics that sets the morphology of crystalline nanoribbons. 
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The electronic properties of crystalline nanoribbons are 
set by their atomic-scale structure, both within the rib- 
bon core and at the edges. The properties are naturally 
modulated by the ribbon morphology and past compu- 
tational studies on graphene nanoribbons (GNRs) have 
identified several shapes for varying widths w and edge 
stresses Tg. For compressive edge stresses, semi-infinite 
sheets exhibit classical edge ripples [TJ |2] that sponta- 
neously twist below a critical width [3]. At even smaller 
widths {w < 1.5 nm), quantum computations show that 
the ribbon becomes flat [4 . Little is known about the 
interplay between geometry and mechanics that drives 
these shape transitions. Intuitively, the accommodation 
of the relative stretch at the edge (or compression for 
Tg > 0) sets the shape stability. More formally, the 
edge stresses induce a non-Euclidean metric as they relax, 
ds'^ = [g{y) dxY + <^y^, where 5 is a surface arc length, x 
and y are coordinates along and normal to the ribbon 
midline and g{ii) is the reference (target) metric that 
captures the interplay; g = I away from the edge and 
changes as we approach the edge in accordance with the 
sign of Tg. The target morphology represents an isomet- 
ric embedding, i.e. no stretching [5 . The standard recipe 
is to seek a shape that minimizes the total energy and 
the metric incompatibility using plate theory which, for 
linear elastic thin sheets (thickness /i <C /, where / is the 
extent of the sheet) and for small strains, reduces to the 
classical Foppl-van Karman (F-vK) equations j6j. 

The atomically thin systems represent the extreme 
limit h ^ and one would naively expect the minimal 
energy surface to be almost inextensible. That is not 
strictly true for the nanoribbons considered here as they 
are inherently confined due to the edge-edge interactions. 
We therefore rely on global minimum energy shapes ob- 
served in high fidelity atomic-scale simulations to identify 
the competing morphologies. The shapes serve as inputs 



to a composite (core-edge) framework, based on the full 
F-vK equations, that is employed to explore the mor- 
phological space as a function of their geometry (w) and 
material parameters (bending and stretching stiffness D 
and S, edge stress and stiffness, Se)- 

The computed shapes are summarized in Fig. [T] In 
addition to the twisted and rippled shapes, we routinely 
observe transversely buckled saddle-like or curled shapes. 
As confirmation, a subset of these shapes have also been 
reported in past studies [1] |3] . A recent study on curled 
ribbons (positive Gaussian curvature) with Tg > shows 
that the morphology is strongly width dependent as the 
ribbon core now deforms out-of-plane to accommodate 
the compressed edges [7 . Interestingly, these shapes can 
also co-exist (Fig. [l^ and Fig. SI). 

In order to quantify the relative stability, we employ a 
composite framework consisting of nanoribbon core, ap- 
proximated as an infinite linear isotropic elastic thin plate 
(width w^ and length / (/i <C i^ <C /)), glued onto ribbon 
edges approximated as elastic strings that stretch or com- 
press in accordance with the edge stress (see Fig. ^). For 
long ribbons, the mid-surface traverse forces are negligi- 
ble and the governing equation for the out-of-plane de- 
flection ( reduces to (see Supplementary Documents) [6] , 
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where T is longitudinal force per unit length (line ten- 
sion). In our framework, an additional incompatibility 
arises at the core-edge interface that relates the ribbon 
force T to the edge stress Tg, 
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where we have set the origin at the midline and r^ are 
the post-buckled edge stresses at the two edges. 






FIG. 1. (color online) (a-b) Shapes observed in atomic- 
scale computations on 10.3nmx2.2nm AGNRs with pristine 
(re < 0) and reconstructed edges (usually Te > 0). (a) Saddle- 
like and twisted shapes for Te < 0. The two shapes also co- 
exist (bottom), (b) Curled shapes for Te > 0. The computa- 
tions were performed using a reactive bond order (AIREBO) 
potential [8 ; see SI for details, (c) Schematic illustration of 
the core-edge composite framework. 



To make contact with these atomically thin systems, it 
becomes necessary to seek solutions in the shallow-shell 
limit. A natural choice for the out-of-plane deflections is 



C{x,y) = -^i^x^ 



Oxy^f{y), 



(4) 



where hz is the longitudinal curvature and 6 is the twist 
per unit length, both assumed constant along the rib- 
bon length. The boundary value problem for the shape 
function f{y) and the ribbon tension T follows from sub- 
stituting Eq. [4] into Eqs. [T] [2] together with free torques 
and shear force balance at the edges. 
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The solution that satisfies Eq. [5] is of the form 

Ci cosh ky cos ky + C2 sinh ky sin ky - 
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Sk (Ci cosh ky cos ky + C2 sinh ky sin ky). 



0. 

(5) 

hC3 

(6) 



and is more conveniently expressed in terms of scaled (di- 
mensionless) variables: curvature R = ^/S/Dw'^hi^ twist 
(per unit length) = ^/S/Dw'^0^ edge stress fg = Tew/D 
and edge modulus, ^e = Se/Sw. In particular, the 
wavenumber kw = \/r/2 and the constant Cs are chosen 

such that the ribbon cannot translate, J^^/2 fiv) dy = 
(Eq. ). The balance between the energy stored in bending 
and stretching the ribbon forms the basis for the post- 
buckled shapes. Setting the energy of a flat relaxed rib- 
bon as our reference, the excess energy 6£ (per unit rib- 
bon length) follows from the shape function and ribbon 
tension (Eq. 85 in SI) and has contributions from bend- 
ing the core S^ = £^ (both mean and Gaussian curva- 
ture) and stretching the core and the edge, £s = ^s +^|- 



Minimizing the ribbon elastic energy with respect to the 
curvature and twist yields the equilibrium shape. 

We first present results for pure twist mainly to bench- 
mark our framework with respect to recent quantum cal- 
culations in AGNRs. Setting /^ ^ in Eq. [6] results in a 
relatively simple expression for the equilibrium twist. 
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It is immediately clear that there exists a critical com- 
pressive edge stress |t*| = 6(1 — i^)(l + 28 e) below which 
the ribbon is flat. The underlying energetics is detailed 
in SI, but the relation can also be reasoned from a simple 
scaling analysis per unit ribbon length. Reverting to un- 
sealed variables, the bending energy scales quadratically 
with the twist angle, S^ = DO^w. The decrease in stretch- 
ing energy arises from the extension of the edge, f| ^ 
—r^O'^w'^^ where rj = re/(l + 2S'e) is the residual edge 
stress following in-plane relaxation. At the onset of buck- 
ling, £b^£s which yields the critical edge stress, fg* ~ 
— {l-\-2Se)' In the limit of large in-plane stiffness Se ^l 
(e.g. GNRs), the critical stress is approximately constant, 
f* ~ —6(1 — i^), i.e. the critical width below which 
the ribbon is flat scales inversely with the edge stress, 
1^* ^ —6(1 — v)D/Te. Past w* the twist angle scales as 
6eq ex (fe — Tg*)^/^, characteristic of a supercritical pitch- 
fork bifurcation. At large widths, the twist decreases 
nonlinearly, Oeq oc (-re/L))^/^^-^/^^-^/^ (Eq. S8). 

Figure [l^ shows the width dependence for pristine AG- 
NRs [v = 0.17, D = 1.5 eV, S = 2000eV/nm2 [1]) with 
compressive edge stresses Tg = — 10.5eV/nm and an edge 
stiffness ^e = 112.6 eV/nm [Hll]. The critical width is 
less than a nanometer, w* ^ 0.8 nm and the twist ex- 
hibits a maximum, Omax ^1.75 deg/A ai w = 1 nm (min- 
imum pitch length Xmin ^ 20.6 nm), which separates the 
two distinct regimes of behavior described above. The 
inset in Fig. [l^ shows the width dependence reported in 
recent quantum computations on AGNRs. The overall 
trends are in excellent agreement with our results, espe- 
cially at small widths, suggesting that, at least in the 
case of AGNRs, the edge-edge interactions are primarily 
elastic in nature. The small deviations (underestimation 
of the critical width and overestimation of the maximum 
twist Omax ^ 1.4 deg/A at w ^ 14 nm) are likely due to 
the differences computational frameworks (classical ver- 
sus quantum). At large widths, the twist decays almost 
linearly in the computations; a similar dependence was 
also seen in separate atomistic computations [3 . This is 
in contrast to the power law behavior and the discrepancy 
is primarily due to the narrow range of widths accessible 
to the computations. 

We now turn our attention to pure bending. The solu- 
tion for the equilibrium curvature hCeq again follows from 
Eq. [6] in the limit ^ ^ 0. The expressions are rather 
lengthy and for the sake of brevity, the width depen- 
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FIG. 2. (color online) Results based on our composite framework, (a) Oeqys.w for pristine AGNRs (re = — 10.5eV/nm). 
(inset) The same variation from quantum computations reported in Ref. 4 and reproduced here, (b) Width dependence of 
longitude curvature Keq for AGNRs with pristine and reconstructed edges, Te = — 10.5eV/nm and Te = 24.5eV/nm respectively, 
(inset) Scaled plot Req/\fe\'^ ys.\re\ in the limit & = for compressive (black) and tensile (red) edge stresses, (c-d) Cross- 
sectional profiles f{y)/wvsy/w for the AGNRs with widths ranging from w = 1 — 10 nm and (c) Te = — 10.5eV/nm and (d) 
Te = 24.5eV/nm. (insets) Schematic illustrations of the ribbon shapes. 



dence is plotted in Fig. ^p for GNRs with varying edge 
stresses. The overall trend is similar to that for Oeq{w) in 
that the flat ribbons spontaneously bend past a critical 
width 1^* with an equilibrium curvature hZeq that rapidly 
increases, exhibits a maximum and then decays non- 
linearly with at large widths. The critical (scaled) edge 
stress for onset of buckling is similar to that for twisting, 
f* = V5{V5u±VQ - u^){1^2Se) (see SI for details) and 
is consistent with simple scaling analysis. For small edge 
deflections d<^w^ the transverse curvature is k^ ~ d/w'^ 
and the curvature preserves the moment balance along 
the ribbon, k, ~ r^S/Dw. Then, the bending energy (per 
unit length) scales as S^ = D{f<i'^ + 5'^ /w^)w. The de- 
crease in energy is due to edge extension/compression 
and scales as E^ = —r^hzS. Equating £b^£e yields the 
scaling, \f^\ ^ (1 + 26'e). The plot shows that the critical 
width for Tg = — 10.5eV/nm is also of the same order as 
that for the twist, to within an angstrom. The quantita- 
tive similarities between the two shapes suggests a small 
energy difference; we address this in detail towards the 
end of this letter. 

The critical curve for Tg = 24.5eV/nm plotted in 
Fig. (2)3 corresponds to reconstructed AGNRs and high- 
lights the effect of edge stress. The equilibrium curva- 
ture in general is larger with a smaller critical width. 
However, the qualitative trends remain unchanged, sug- 
gesting a more general interplay between hZeq and Tg dur- 
ing post-buckling. As confirmation, the inset in Fig. (2)3 
shows the scaled plot, Req/\^e\'^ vs. |fe|, predicted by our 
framework in the limit Se = 0. Positive edge stresses 
result in curled ribbons with larger critical width - the 
difference is entirely due to the Poisson's ratio (see SI). 
At large widths, Req ex \fe\^ with the scaled exponent 
n = 4/3. The unsealed curvature therefore exhibits a 
power-law decay with width, tv^q ex k;~^/^, consistent 
with recently reported results for GNRs with Tg > 0. 
At intermediate widths in the vicinity of the peak cur- 



vature, the exponent oscillates about n = 4/3 (see SI). 
The deviations in the scaling are associated with tran- 
sitions in the ribbon profile. Figure [Sj^-d shows these 
profiles for saddle-like and curled AGNR profiles as a 
function of width. Note the excellent agreement be- 
tween our framework and computations. The profile for 
w = 9.9 nm is representative of the shapes at large widths 
in that the out-of-plane deflection are strongly localized 
at edges that do not interact. The edge boundary layer 
is preferentially strained and therefore bent, resulting in 
a double- well profile with a relatively flat midsection. At 
small widths {w = 2.2 nm and 4.6 nm), the deflection 
penetrates through the width due to the edge-edge inter- 
actions mediated by a strained core as well as the bound- 
ary conditions. The profiles evolve such that the mean 
curvature hz -\- hz^ decreases with width. 

Compressive edge stresses also lead to co-existing 
twisted and saddle-like shapes (Fig. Ilk) . Insight into the 
relative stability can be gleaned from the athermal be- 
havior, readily available through our framework for any 
combination of 6 and hz (see SI). The results are plotted in 
Fig. |3^ as a contour plot of the excess energy S£ for AG- 
NRs with Tg = — 10.5eV/nm and fixed width, w = 3nm. 
The colored regions represent energy gain associated with 
the morphology, i.e. S£ < with respect to planar yet 
relaxed ribbons. The origin {0 = K: = 0) is an unsta- 
ble local maximum, as expected. The local minima for 
pure twist and bending (solid circles) correspond to the 
equilibrium points in Figs. [2^ and[2]3. Pure bending is fa- 
vored, yet the energetic cost associated with a co-existing 
twist is exceedingly small (0.01 eV/nm) and can be easily 
overcome by the thermal energy available at room tem- 
peratures {ksT ^ 0.025 eV/atom). 

The width dependence of the energy gain associated 
with both morphological classes yields a comprehensive 
picture of the shape transitions. Figure ^jp shows such a 
plot for GNRs with Tg = ±10.5eV/nm. The individual 
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FIG. 3. (left) Contour plot of the elastic energy change 68 predicted by our framework for bent and twisted AGNRs with 
pristine edges (re = — 10.5eV/nm). (middle) Width dependence ^^(w) for minimal shapes of pristine and reconstructed AGNRs 
(re = itl0.5eV/nm). The energy for edge rippling (solid black line) is based on wavelength observed in simulations and an 
assumed form for the shape 1 . See text and SI for details, (c) (right, top) The minimal shape for an 18.5 nm long tapering 
ribbon with opening angle is a = 7r/6. The vicinal edges are at an angle 7r/12 off the armchair orientation. The color indicates 
the longitudinal curvature k,. (right, bottom) Variation in the ribbon curvature k,{s) towards the tapering end. 



contributions are plotted in detailed in SI. For compres- 
sive edge stress, the two shapes are energetically similar 
for w* < w < 3nm; the behavior is consistent with our 
observations so far. Pure bending becomes more favor- 
able beyond a transition width w = 3nm. In the large 
width limit, the energy gap that stabilizes the saddle- 
like morphologies scales as 5£{0 = 0) — S£{hi = 0)~ 
T^/{Sw){\fe\'^^^ — I). The non-linear dependence on the 
edge stress offers insight into the relative shape stability 
in general; small compressive edge stresses (e.g. hydro- 
genated GNRs) stabilize the twist such that the transi- 
tion to saddle shapes shifts to larger ribbon widths, while 
large stresses favor saddle shapes. 

At large widths, the edges will eventually ripple [10]. 
While a general treatment of this transition is beyond the 
scope of the present study, the transition width for the 
specific case of AGNRs can be estimated based on the 
equilibrium wavelengths observed in prior computations 
(A = 5.6 nm) (TJ |3] and an assumed shape of the ripples 
(see SI). Note that the corrections due to the edge-edge 
interactions are absent in our estimate. Based on the ap- 
proximate excess energy of a ribbon so rippled (Fig. Isb), 
we expect a twist-rippling transition dit w = 6 nm and a 
bending-rippling transition at widths of the order of tens 
of nanometer. The transition width should again scale in- 
versely with the edge stress. For example, computations 
on AGNRs with Tg ~ — 26eV/nm show that edge ripples 
become favorable relative to twist (the saddle shapes were 
ignored) at widths as low as 1.5 nm [3 . Finally, in the 
case of tensile edge stresses, the twisted morphology is 
obviously not viable. Interestingly, the energy associated 
with curled shapes is higher than that for saddle-shapes 
for the same width. At much larger widths, the com- 
pressed midline buckles into dimples that decay towards 
the edges, analogous to the bending-rippling transition. 

In conclusion, the interplay between ribbon width and 
mechanics that sets the shape highlights the importance 
of structural and chemical control of the edges. It also 
raises intriguing questions on shapes in ribbons with local 



variations in geometry and/or edge properties. As a case 
in point, Fig.jS]^ shows the computed minimal shape of a 
tapering GNR [11 . The shape changes towards the ta- 
pering end in accordance with our predictions; the wider 
end is rippled, settles into a saddle-like curved shape, and 
then flattens out into a slightly twisted tapering end. The 
shape is not unlike the curled tips of growing leaf blades 
and flower petals, highlighting a general principle where 
geometrical (and possibly material) modulations are nat- 
urally amplified as controllable shape transitions. 
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